Combining the advantages of superconducting MgB 2 and CaC@ in one material: 
suggestions from first-principles calculations 
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We show that a recently predicted layered phase of lithium monoboride, Li2B2, combines the key 
mechanism for strong electron- phonon coupling in MgB2 (i.e., interaction of covalent B a bands 
with B bond-stretching modes) with the dominant coupling mechanism in CaC6 (i.e., interaction 
of free-electron- like interlayer states with soft intercalant modes). Yet, surprisingly, the electron- 
phonon coupling in Li2B2 is calculated to be weaker than in either MgB2 or CaC6- We demonstrate 
that this is due to the accidental absence of B ir states at the Fermi level in Li2B2. In MgB2, the ir 
electrons play an indirect but important role in strengthening the coupling of a electrons. Doping 
Li2B2 to restore ir electrons at the Fermi level is expected to lead to a new superconductor that 
could surpass MgB2 in T c . 

PACS numbers: 74.25.Jb, 74.25. Kc, 74.70.Ad 



I. INTRODUCTION 

The discovery of superconductivity in MgB 2 at T c — 39 
K led to a new paradigm for achieving high critical 
temperatures in conventional superconductors, namely 
strong electron-phonon coupling due to metallization of 
strong covalent bonds. 1 In MgB2, this is achieved through 
the interaction of electrons forming covalent B-B bonds 
in the graphene-like layers with vibrational modes that 
stretch and compress these bonds. 2 . Unfortunately, ef- 
forts to find related materials with even higher T c have 
not been fruitful. This is basically because the electronic 
states that form the covalent bonds in MgB2 are quasi 
two-dimensional, which makes it difficult to enhance their 
density of states and their coupling to phonons through 
simple doping. 3 Significant progress may only be possible 
through the introduction of different types of electronic 
states at the Fermi levels 

The recent discovery of superconductivity in CaC6 at 
T c = 11.5 Kj^ the highest critical temperature among 
graphite intercalation compounds (GICs), offers an al- 
ternative path to increasing T c . Although CaC6 is struc- 
turally similar to MgB2, it appears that the key elec- 
trons and phonons responsible for its superconductivity 
are distinct from those in MgB2^^ In particular, free- 
electron-like states that are concentrated in the region 
between graphene layers dominate the electron-phonon 
interaction through coupling with soft intercalant modes 
and with bond-bending modes of the graphene sheet. In- 
cipient lattice instabilities associated with the soft inter- 
calant phonon modes, however, may present obstacles to 
attaining significantly higher T c in GICs. 5,7 

It is intriguing to ask whether the MgB2 and CaCg 
mechanisms for strong electron-phonon coupling can be 
incorporated into a single compound to create a super- 
conductor with the potential to surpass MgB2. Although 
neither MgB2 nor CaCg can be easily modified so as to 
incorporate the key electronic and phonon players from 
the other material, we show here that it is indeed possible 



to combine the MgB2 and CaC6 paradigms in a different 
material. Such a material has a good chance of exceeding 
both MgB 2 and CaC6 in terms of T c . 

The material of interest is a hypothetical layered form 
of lithium monoboride, Li2B2, that recent calculations 
have predicted to be competitive in energy with known 
LiB phases.^ The B graphene sheets in Li2B2 are nearly 
identical to those in MgB2. At the Fermi level, there are 
a bands derived from B p xy orbitals that form covalent 
bonds within the layers, as in MgB2, and there are inter- 
layer free-electron-like bands, which we shall call £ bands, 
as in CaC6- By accident, the B p z -derived ir states that 
are present in both MgB2 and CaCg are missing at the 
Fermi level in Li2B2. Based on the similarity of the in- 
plane B physics with MgB2, Kolmogorov and Curtarolo 
have suggested that this material could be a good super- 
conductor in the mold of MgB2^ 

In this paper, we analyze the electronic structure of 
Li2B2 and show that the electron-phonon interaction has 
sizable contributions not only from the a bands that 
dominate the coupling in MgB2, but also from the inter- 
layer £ bands, which are the key players in CaCg. Sur- 
prisingly, the total electron-phonon coupling in Li2B2 is 
found to be weaker than in either of the other two ma- 
terials. We show that this is precisely because of the 
missing ir electrons at the Fermi level. Though the 7r 
bands in MgB2 and CaCg give only moderate contribu- 
tions to the electron-phonon coupling, 7r electrons can 
indirectly strengthen the coupling of a electrons with B 
bond-stretching phonon modes through screening. Dop- 
ing Li2B2 to add ir sheets to the Fermi surface would 
enhance the a contribution to the electron-phonon cou- 
pling, as well as add a modest ir contribution. While a 
relatively low T c is predicted for pristine Li2B2, an ap- 
propriately doped material could have a T c higher than 
even MgB2. 
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II. CALCULATIONAL METHODS 

Our analysis is based on density functional calcula- 
tions. The electronic structure was calculated using the 
WIEN2K general potential LAPW code,- 10 with the GGA 
functional in the PBE fornix Local orbitals were added 
to relax the remaining linearization error and APW or- 
bitals were used to improve the convergence with re- 
spect to the cutoff parameter RK ma ^. In most calcula- 
tions RKm&x — 7 was used, except for the convergence 
tests with RK mSiX = 8. Up to 918 irreducible k-points 
were used (32x36 x 16). The phonon spectrum and 
electron-phonon coupling functions were calculated using 
the linear-response method 12 within the local density ap- 
proximation. Norm-conserving pseudopotentials^ were 
used, with a plane-wave cutoff of 50 Ry. The dynami- 
cal matrix was calculated on a grid of 6 x 6 x 2 phonon 
wavevectors q in the Brillouin zone and Fourier inter- 
polated to denser meshes for calculation of the phonon 
density of states (DOS). For calculation of the electron- 
phonon matrix elements, the electronic states were sam- 
pled on grids of 36 x 36 x 12 k-points in the Brillouin 
zone. 



III. RESULTS AND DISCUSSION 

The structure of the hypothetical layered compound 
Li2B2, first proposed in Ref. [S| (in which it is called 
MS2), consists of honeycomb B sheets, each surrounded 
above and below by layers of Li on a triangular lattice. 
There is a horizontal shift between neighboring Li-B-Li 
sandwich units, leading to an ABAB... stacking sequence 
of the Li-B-Li trilayer units. While we find that the opti- 
mized structural parameters, particularly the interlayer 
distances, are sensitive to the details of the electron-ion 
potential and to the exchange-correlation functional, the 
band structures calculated within the different approxi- 
mations are very similar. For the results presented here, 
we have fixed the structural parameters to the optimized 
values reported in Ref. H. 

The band structure of Li2B2 near the Fermi level is 
shown in Fig. [T] As expected, the a (blue) and ir (green) 
states are derived primarily from in-plane B p x ^ y and out- 
of-plane B p z orbitals, respectively. The C bands (ma- 
genta) have substantial Li sp and B p z character. The 
charge density of the £ electrons, plotted in Fig[2l is con- 
centrated in the region between Li-B-Li trilayer units, 
though hybridization with B p z orbitals is evident. This 
is similar to the character of ( bands in GICs£ i 14 ' 15 The 
Fermi level in Li2B2 lies at the crossing of the ir and ir* 
bands at K, reminiscent of pure graphite. Note, however, 
that in graphite, the coincidence of the Fermi level with 
the 7r-7r* crossing is mandated by the absence of other 
bands near Ep, while in Li2B2, it is accidental. (Expand- 
ing c/a by ~ 10%, for example, moves the tt-tt* crossing 
away from the Fermi level by more then 0.5 eV.) Com- 
pared to pure graphite, the presence of Li atoms between 



4 
3 
2 
1 

> 

LLI 

-1 

-2 
-3 

-4 



FIG. 1: (Color online) Li2B2 band structure near the Fermi 
level {Ef = 0). The a, n, and £ bands are plotted with blue, 
green, and magenta symbols, respectively. 
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FIG. 2: (Color online) Isosurface of Li2B2 charge density cor- 
responding to the £ state near the Fermi level, half way be- 
tween the r and A points. The top, center, and bottom plane 
in the figure contain honeycomb sheets of B. The charge den- 
sity in the center B plane is shown. Two Li layers lie between 
each pair of B sheets. The £ states are primarily concentrated 
in the interlayer region, but have some B p z character as well. 



the B-based graphene sheets i) makes the ir bands more 
3D and lowers their energy relative to the a bands, and 
ii) lowers the C bands due to the attractive ionic potential 
in the interlayer region where the electronic charge in the 
£ bands is concentrated*!^ 4 - With these relative shifts in 
the bands, the a bands, which are completely filled in 
pure graphite, become partially occupied in Li2B2, while 
the £ bands, which are unoccupied in graphite, cross the 
Fermi level in Li2B2. Hence Li2B2 can be viewed as a 
self-doped analog to pure graphite, with holes in the a 
bands exactly compensated by electrons in the £ band. 
In contrast, in CaCg the it* band is partially occupied as 
well. 

In comparing the band structures of Li2B2 and MgB2, 
which both have B honeycomb sheets separated by "in- 
tercalant" layers that donate one electron per B site, key 
differences in the relative energies of the different types 
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FIG. 3: (Color online) Fermi surface of Li2B2. The cylindrical 
sheets centered on the F-to-A line arise from the 2D a bands, 
while the ellipsoids are due to free-electron-like £ bands. 



of bands can be attributed to the difference in the sep- 
aration between B layers. With two Li layers between 
B sheets in Li2B2, the B-B interlayer distance in Li2B2 
(5.4 A) is significantly larger than in MgB2 (3.5 A). This 
makes the p^-derived a bands more 2D in Li2B2, but 
has little effect on their in-plane dispersion. As in MgB2, 
these bands give rise to nearly cylindrical hole sheets of 
the Fermi surface surrounding the V to A line. As can be 
seen in the Fermi surface plotted in Fig. [3j there are two 
nearly degenerate pairs of a cylinders. The larger sepa- 
ration between B planes also causes the B p z -derived tt 
bands to be more 2D in Li2B2 and raises them in energy 
so that the tt* band that is partially occupied in MgB2 be- 
comes unoccupied (and just touches Ef at certain points 
in the Brillouin zone). The Fermi surface due to the tt 
bands thus disappears in L12B2. On the other hand, the 
free-electron-like £ bands, which are unoccupied in MgB2, 
drop in energy when the interlayer distance is increased 
since there is a larger interlayer volume available for the 
electrons to occupy In Li 2 B 2 , the £ bands produce a 
slightly squashed ellipsoidal Fermi surface around T, with 
a smaller effective mass for motion along the z direction 
than in the plane. These ellipsoids are huge compared to 
the small ( balls in the CaCg Fermi surface^ even though 
CaC6 has the most occupied £ band among all GICs. 14 

The topological simplicity of the Fermi surface of 
MgB2, in which sheets of the Fermi surface correspond- 
ing to bands of different character are widely separated 
in reciprocal space, is lost in Li2B2. Not only do the £ 
ellipsoids cross all four a cylinders, they also cross each 
other. At all of these crossings, there is only weak hy- 
bridization between the bands. The multiple crossings 
between sheets of the Fermi surface arising from bands of 
different character, along with the presence of bands that 
nearly touch the Fermi level, make Brillouin-zone inte- 
grations sensitive to k-point sampling. This is especially 
an issue for calculations of the electron-phonon coupling 
constants, which involve double integrals over the Fermi 
surface. The quantitative electron-phonon coupling re- 
sults presented here likely have larger uncertainties than 
usual but this does not affect the qualitative picture that 
emerges. 

The phonon density of states F(lS) and electron- 
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FIG. 4: Li2B2 phonon density of states F{ui), electron-phonon 
spectral function a 2 F(uj), and ratio of a 2 F to frequency lo. 
The electron-phonon coupling constant A, which is propor- 
tional to the inverse-frequency-moment of a 2 F, has significant 
contributions from both low-frequency and high-frequency 
modes. 



phonon spectral function a 2 F(ui) calculated for Li2B2 
are shown in Fig. 21 The Ei g in-plane bond-stretching 
B modes that are the key players in superconductivity of 
MgB2 are also important in Li2B2, as indicated by the 
peak in the spectral function near 80 meV. In both mate- 
rials, these vibrational modes couple strongly to the B a 
bands. Given the similarity of the B in-plane physics in 
the two materials, one expects similar deformation poten- 
tials for the Ezg modest Indeed, our frozen-phonon cal- 
culations using LAPW-GGA produce exactly the same 
deformation potential at T : a displacement u opens a 
gap As such that As/u = 26 eV/A, precisely the same 
as found for MgB2.— Furthermore, while it is difficult to 
separate the electronic DOS in Li2B2 into contributions 
from different bands, as is routinely done in MgB2, the ef- 
fective masses derived from the dispersion of the a bands 
are within 10% of those in MgB2 (the heavy band is 9% 
heavier and the light band 8% lighter). If the phonon 
frequencies were the same, one would then expect the 
net intraband electron-phonon coupling constant for the 
a bands, A CTCr , to be practically the same as in MgB2. 
Surprisingly, the phonon frequencies differ considerably. 
At the zone center, the Ei g phonon frequency of 81 meV 
in Li2B2 is signficantly higher that the corresponding 67 
meV in MgB2r' 1 ^ which likely reduces \ aa from ps 1 to 
ps 0.6 -0.7. 

The hardening of the in-plane B bond-stretching 
modes in is related to the loss of tt electrons at the 

Fermi level. The same electron-phonon coupling process 
that gives rise to superconductivity also leads to elec- 
tronic screening of phonons. For a zone-center phonon 
mode, the softening due to screening by metallic elec- 
trons is given by 



Au? = -4uj(g 2 )N(0), 



(1) 



where (g 2 ) is a Fermi-surface average of the square of 
the electron-phonon matrix element g, and N(0) is the 
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electronic DOS at the Fermi level^^ Because the ma- 
trix clement g contains a factor of l/\/uJ, the right-hand 
side of Eq. Q] is independent of the phonon frequency 
to. In MgB2, both the a and tt electrons participate in 
screening the bare E^ g phonons. For the Ei g mode in 
MgB2, it has been estimated that cu 2 is reduced by about 
1500 meV 2 due to screening by tt electrons, and by about 
2000 meV 2 due to screening by a electrons, and that the 
unscreened frequency is luq ~ 90 meV. 2 In Li 2 B 2 , the 
unscreened frequency and the softening due to a bands 
should be roughly the same magnitude as in MgB 2 , which 
gives a rough estimate for the screened Ei g phonon fre- 
quency of u> = \Jijj 2 — (Acj 2 ) ct w 78 meV. This is consis- 
tent with our actual calculations for Li 2 B 2 , which yield 
a frequency of 81 meV at the zone center. Since only 
phonons with q xy < 2k F , where k F is the radius of the 
cylindrical a sheets of the Fermi surface, can couple to 
<j electrons, softening by a electrons is only significant 
when the in-plane component of the phonon wavevec- 
tor is less than 2k p. Indeed, at the zone boundary, our 
calculations yield a frequency of 92 meV for the lowest 
E2 g branch, consistent with estimates for the unscreened 
E 2g frequency. While Li 2 B 2 also has £ electrons, which 
are not present in MgB 2 , these are largely localized in 
the interlayer region and do not couple strongly with the 
B bond-stretching modes, so they are less effective in 
screening the Ei g phonons. The harder Ei g frequencies 
in Li 2 B 2 compared to MgB 2 can thus be attributed to 
the absence of screening by tt electrons. 

In addition to the high-frequency Ei g peak in the 
spectral function, there is a low-frequency peak in a 2 F 
arising from ultrasoft modes that involve rigid shifts of 
the B planes in the z direction coupled with both in- 
plane and out-of-plane motion of the Li atoms. The 
frequency of these ultrasoft modes is sensitive to cal- 
culational details, so at T, for example, the harmonic 
frequency of about 8 meV obtained in our LDA linear- 
response calculation is somewhat lower than the esti- 
mate of 11 meV based on force-constant methods us- 
ing PAW GGA calculations. 8 In any case, these low- 
frequency normal modes couple with the C bands, sim- 
ilar to the situation in CaC6, where Ca sliding modes 
interact strongly with interlayer electrons. Although rel- 
atively small, the low-frequency peak in a 2 F gives a sim- 
ilar contribution to the isotropic electron-phonon cou- 
pling constant A = 2 J duia 2 F/uj as the high-frequency 
peak (see Fig. |U bottom panel). The total A is calcu- 
lated to be 0.57, with the low-frequency peak contribut- 
ing about 0.15, and the Ei g peak contributing about 
0.2. The logarithmically averaged phonon frequency of 
u)\ n = exp[A _1 J duohi{uj)a 2 F(uj)/u)] = 39 meV, is well 
below the average frequency of uj ave = 59 meV, reflect- 
ing the importance of the ultrasoft modes. 

From the point of view of T c , the prospects for Li 2 B 2 do 
not look optimistic compared to MgB 2 , because of both 
the lack of a 7r-band contribution to A, and, perhaps more 
importantly, the hardening of the key Ei g phonons that 
are no longer screened by tt electrons. Neglecting possi- 



ble multiband effects and using the McMillan equation^ 
with a typical Coulomb parameter of fj,* = 0.12, we es- 
timate T c « 7 K. In comparison, the estimated isotropic 
T c in MgB 2 is 25-30 K. Enhancement of T c due to inter- 
band anisotropy is expected to be weaker in Li 2 B 2 since 
the two groups of electrons contribute more equally to 
the coupling with phonons. But even if it were of the 
same order as in MgB 2 , that is, ~ 30% in the effective A, 
it would only increase T c to about 15 K. The bottom line 
is that Li 2 B 2 , if synthesized, might be competitive with 
CaCg in terms of T c , but not with MgB 2 . 

However, that is not the end of the story. 

While the undoped Li 2 B 2 is isoelectronic with MgB 2 , 
and doping of MgB 2 is known to unequivocally depress 
T c , the same need not be true for Li 2 B 2 . Electron dop- 
ing of MgB 2 lowers the a DOS, and hence the electron- 
phonon coupling, since the bands are not ideally 2D. 
(Though attempts have been made, hole doping of MgB 2 
has not yet been achieved.) In Li 2 B 2 , the a bands are 
much more 2D, so doping should have little effect on their 
DOS and electron-phonon coupling. On the other hand, 
electron doping would start filling the tt bands, increas- 
ing their DOS. Assuming rigid bands, doping on either B 
or Li sites could be used to raise the Fermi level by nearly 
1 eV without losing the contribution from the a bands. 
We estimate that this would increase the 7r-band DOS to 
roughly half the value in MgB 2 . Doping on the Li sites 
(with Al, Be, Mg, etc.) should yield even higher tt DOS, 
since the extra ionic charge in the Li planes would lower 
the tt bands relative to the a bands. Experience shows 
that it is hard to dope more holes into the a bands than 
there already are in pure MgB 2 . Since Li 2 B 2 starts with 
more holes in the a bands than MgB 2 , it might even be 
easier to synthesize electron-doped Li 2 B 2 than the un- 
doped material. Assuming that the tt bands in doped 
Li 2 B 2 couple with the Ei g modes at about the same level 
as in MgB 2 , this additional coupling will soften the Ei g 
phonons and enhance their coupling with the a bands, 
without strongly affecting the coupling of £ bands with 
soft modes. This is an extremely favorable combination 
that is likely to have a higher T c than MgB 2 . To a lesser 
extent the same holds if the accidental alignment of the 
Fermi level with the tt-tt* crossing is modified by uniaxial 
pressure. 

While Li 2 B 2 and related layered lithium borides pro- 
posed by Kolmogorov and CurtaroloS are only hypo- 
thetical materials as yet, their calculations of the Li-B 
phase diagram have suggested a route for synthesizing 
these materials using pressured We have shown that al- 
though Li 2 B 2 seems to combine the key mechanism for 
strong electron-phonon coupling in MgB 2 with the domi- 
nant coupling mechanism in CaCg, each of which has the 
highest superconducting transition temperature in its re- 
spective class of materials, the transition temperature of 
Li 2 B 2 is expected to be lower because of the lack of the 
seemingly unimportant tt electrons at the Fermi level. An 
appropriately doped version of the material, in which tt 
electrons are restored at the Fermi level, however, could 
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result in a promising new superconductor that combines 
the paradigms of MgB 2 and CaC@ in a single material 
without loss. Further calculations using the virtual crys- 
tal approximation or the supercell approximation to ex- 
plore the quantitative effects of doping on the electron- 
phonon coupling in this material would be of interest. 
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